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Abstract 

The significance of the recent finding of the velocity distribution function of 
the steady-state Boltzmann equation under a steady heat current obtained 
by Kim and Haykawa (J. Phys. Soc. Jpn. 72, 1904 (2003)) is discussed. 
Through the stabihty analysis, it seems that the steady solution is stable. 
One of possible applications to the nonequiliburium Knudsen effect in which 
one cell at equilibrium is connected to another cell under the steady heat 
conduction is discussed. This solution apparently shows that steady-state 
thermodynamics proposed by Sasa and Tasaki cannot be used in a naive 
setup. The preliminary result of our simulation based on molecular-dynamics 
for nonequilibrium Knudsen effect is also presented to verify the theoretical 
argument. 

I. INTRODUCTION 

Kinetic theory of dilute gases has played crucial roles in the history of nonequilibrium 
statistical mechanics. The theoretical treatment of linear transport problems is well estab- 
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lished including the derivation of Navier-Stokes equation and the relation between Boltz- 
mann equation and the linear-response theory through a number of studies of Boltzmann 
equation [1-3]. 

The transport problems under the nonlinear-nonequilibrium circumstance, however, is 
still controversial. In fact, we know Burnett's classical work in which he obtained an approxi- 
mate hydrodynamic equation at the second order nonequilibrium condition [4] , but we realize 
that hydrodynamic equation he obtained is unstable for perturbations [5,6]. Therefore, one 
may suggest that the validity of Chpaman-Enskog method is so limited that we should use 
Grad expansion [7] to discuss the higher order perturbations such as Burnett order and 
super-Burnett order [6]. On the other hand, when we obtain the steady solution of Burnett 
equation, the stability of the steady inhomogeneous solution may be different from that 
of dynamical Burnett equation for homogeneous perturbation. Therefore we cannot judge 
whether the calculation at Burnett order based on Chapman-Enskog method is useless even 
when Burnett's hydrodynamic equation is unstable. Unfortunately, these difficulties might 
be regarded as unimportant, because the solution at Burnett order is practically useless 
except for few cases such as the analysis of shock waves. 

In these days, however, the solution of Boltzmann equation at Burnett order has at- 
tracted much interest among researchers who try to extend the linear nonequilibrium sta- 
tistical mechanics to nonlinear nonequilibrium statistical mechanics [8,9]. In fact, since all 
of attempts are formulated to match them with the linear response theory, their consistency 
with nonlinear theory like Burnett solution should be important. 

Recently, Kim and Hayakawa [10] have extended Burnett's classical work on Boltzmann 
equation and obtained the explicit steady solution of Boltzmann equation under the steady 
heat conduction at second (Burnett) order. The accuracy and numerical stability of this 
solution is confirmed by Fushiki [11] from his molecular-dynamics simulation. This stable 
behavior does not contradict with blown-up tendency of Burnett equation for perturbation 
of short wavelength or high frequency, because the simulation does not involve any pertur- 
bation with short wavelength or high frequency. Their solution is useful [12] to examine the 
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validity of nonequilibrium theories such as information theory [8] and the steady-state ther- 
modynamics (SST) [9] based on thermodynamic arguments. They also succeed to apply the 
solution to calculate the rate of chemical reaction under the nonequilibrium circumstance in 
which the linear nonequilibrium effects are canceled from the consideration of symmetries 
[13]. 

The purpose of this paper is to (i) summarize the achievement of Kim and Hayakawa, 
(ii) to discuss the stability of the steady solution at Burnett order, and (iii) to introduce 
the result of the molecular-dynamics (MD) simulation for nonequilibrium Knudsen effect to 
examine SST. 

The organization of this paper is as follows. In the next section, we briefly summarize 
the outline of the solution obtained by Kim and Hayakawa [10]. In section III, we will 
discuss the stability of the solution. In section IV, we will discuss Knudsen effect under a 
nonequilibrium condition based on MD simulation. 

II. THE OUTLINE OF THE SOLUTION BY KIM AND HAYAKAWA 

In this section, we briefly summarize the solution of steady Boltzmann equation under 
the heat conduction obtained by Kim and Hayakawa [10]. Suppose a system of dilute gas 
in a steady state whose velocity distribution function is /i = /(r,vi). The steady state 
Boltzmann equation is written as 

Vl-V/l = J(/l,/2), (1) 

where the collision integral J(/i, f2) for hard core molecules is expressed as 

Afi, h) = J I J if[f2 - fif2)gbdbdedv2. (2) 

Here the velocity distribution /i and /2 change to /{ and /2 by a binary collision. The 
relative velocity of the two molecules before and after the interaction has the same magnitude 
g — |vi — V2I. The relative position of the two molecules is represented by the impact 
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parameter b and e which represents the angle of the plane to which g and g'. The impact 
parameter b is identical to explicitly determined by specifying the kind of molecules. 
The velocity distribution function /i can be expanded as: 

A = /f ^ + fi'' + + ■■■ = f?\i + + <f>? + ■■■), (3) 

where the small parameter of the perturbation is Knudsen number , i.e. mean free path 
of the molecules is much less than the characteristic length of the change of macroscopic 
variables, fi*^ in eq.(3) is the local Maxwellian distribution function written as 

where m is mass of the molecules. n(r) and T(r) denote the density and the temperature at 
position r, respectively. Substituting eq.(3) into the steady state Boltzmann equation (1), 
we arrive at the following set of equations: 

L[A'M'-^.-yfr\ (5) 

for 0{K) and 

Lm<t>f'-^.-Vf?-J{f?fi\ (6) 
for 0{K'^)^ where the linear integral operator L[ff'^] is defined as 

L[/f ^ ( f ( fi'^f^'^ {X[ -X, + X',- X,)gbdbdedv,. (7) 



It is important to consider the solubility conditions of these integral equations, for the 
steady state Boltzmann equation, which are given by 

/ ■ v/f ^dvi = 0, (8) 

for 0{K), where $i is one of the coUision invariants: 

$1 = 1, $2 = "iVi, $3 = ^mv^. (9) 



Substituting eq.(4) into the solubility condition (8), it can be shown that nT is uniform. We 
use this result in our calculation. The solubility condition for 0{K^) is given by 

/ ^^vi • Wfi'Uvi = 0. (10) 

To make the solutions of the integral equations (6) and (7) definite, five further conditions 
must be specified; we identify the density: 

n(r) = J hdw, = j fi'Uv,, (11) 

the temperature: 

5!!Mn£)=|!!^;.,v. = |!!^;<«W, (12) 

and the mean flow: 

n(r)u(r) = J v^f^dv, = J v^f^'^d^r,. (13) 

Here we assume that no mean flow u = exists in the system. 

Skipping the derivation of the solution at Burnett order, let us write the explicit form of 
the steady solution for hard core molecules as under the temperature gradient: 

+ Y. '■!62,(24 - 4 - cl)r(r + |)S|(c2)}], (14) 

r>0 ^ 

where 61,., 6or a-nd for r < 7 can be determined as listed in Tables.I. Note that we have 
confirmed the convergence of series b^r by changing r. Here represents x— component of 
the heat-flux in the cell on the right hand side in the nonequilibrium steady state which 
corresponds to the x component of the heat flux in eq.(17). In eq.(14), r{x) is the Gamma 

function, c = \^Jm/2T, and Sl{X) is the Sonine polynomial defined by 

00 

(1 - u;)-'-'e-^ = J2T{p + k + l)Sl{X)u^. (15) 

p=0 
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This kind of analysis, of course, is possible to two-dimensional dilute hard-core gases. Kim 
[14] has already obtained the result to compare his result with MD simulation. 

In the derivation the solubility conditions for 0{K^) is important. Thus, eq.(lO) leads 
to the condition: 

V • q(^) = 0, (16) 
where q^^\ i.e. the heat flux for f[^^ can be obtained as 

with bn listed in Table.I. It must be emphasized that q^^^ should satisfy 

q(') = (18) 

from the symmetry consideration of /(^). Thus, the solubility conditions for /^^^ of the steady 
state Boltzmann equation lead to the fact that the heat flux is constant to the second order. 
From eqs.(16) and (17), we also obtain an important relation between (VT)^ and V^T as 

^ + = 0. (19) 

From eq.(19), the term of V^T can be replaced by the terms of (VT)^. 

Amongst several applications of the solution, in this paper, we emphasize that the so- 
lution can examine the vahdity of SST proposed by Sasa and Tasaki [9]. They propose a 
non-trivial nonequilibrium effect in a simple setup as follows( Fig. 1). There are two cells 
whose size along x axis is Lq. Both of the cells are filled with dilute gases and connected to 
each other by a small hole at x — whose linear dimension along x axis and its diameter 
are less than the mean free path. The cell on the left hand side is at equilibrium of tem- 
perature To, while the cell on the right hand side is in a nonequilibrium steady state under 
a temperature gradient caused by the right wall at temperature Ti at x — Lq and the thin 
mid- wall at temperature Tq. 
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FIGURES 




x=0 



FIG. 1. Setup of nonequilibrium Knudsen effect. 



Sasa and Tasaki predicted the following results [9]. When the system shown in Fig.l is 
in a steady state, the osmosis AP 



= Pxx ~ Pqi 



(20) 



always positive, where Pq is the pressure in the equilibrium cell and Pxx is the xx-component 
of the pressure tensor in the nonequilibrium cell which is defined by 



Pij = j dvm{vi - Ui){vj - Uj)f{v). 



(21) 



In addition, there is a relation among n{0), P^x, no and Pq as 

n(0) fdP,x\ 



(22) 



^0 ydPojT^^j 

where n{0) is the density of the cell in the nonequilibrium steady state around the wall and 
no is the density of the cell at equilibrium. 

A nonequilibrium steady state is achieved when the mean mass flux at the hole is zero. 
Neglecting the Knudsen layer effect, i.e. the shp effect around the 'wall' [3,15], we can 
calculate the balance of mass flux as 



/ dvx / (fv±mVxfo + dvx d^v±mVxf\x=o = 0, 

Jo J J-oo J 



(23) 
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where v± represents the components of the velocity which are orthogonal to Vx, i.e. Vy 
and Vz- Prom eq.(23), we can obtain the relation between the density of the cell in the 
nonequilibrium steady state around the hole n(0) and that of the cell at equilibrium no as 

n(0)=no[l + A„^], (24) 

to the second order, where A„ ~ 8.21 x 10~^ in the approximation until 7th order expansion 

of bkr- 

Similarly we can calculate pressure tensor as 

■ ■ IT) 

P^J ^ nT[S,j + X'^-^J, (25) 

with the unit tensor dij and the number tensor A|f = —4.600 x 10~^ and Xp = 2.300 x 10~^. 

Prom substitution of eqs.(24) and (25) into eq.(22) leads to the relation Xj?/Xn — —2 
in SST. On the other hand, our results conflict with this relation: our results become 
Af?/A„ = —0.5604 for hard core molecules. This is a negative result to SST, though the 
treatment of the wall in our calculation is too primitive. We also indicate that our treatment 
in eq.(24) assumes that we can use the local distribution function because of small Knudsen 
number, though, for dilute gases with finite temperature gradient, the place with the mean- 
free path away from the hole has different temperature from that at the hole because Knudsen 
number should be finite. 



III. STABILITY OF THE SOLUTION 

In this section, let us discuss the stability of the second order solution in eq.(14). Honestly 
speaking, it is difficult to obtain the complete stabihty analysis of the problem. Therefore, 
here, we focus on the stability of five hydrodynamic modes which are zero-eigen modes of 
linearized Boltzmann operator and most dangerous ones in the stability. 

We should recall that the steady solution with the temperature field T{x) is obtained 
under the assumption of the zero mean (hydrodynamic) flow. With the aid of eq.(18) one of 
the hydrodynamic equations corresponding to the solution is the steady diffusion equation: 
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V.(.VTW)=0, . = (26) 

where k is the heat conductivity. Another hydrodynamic equation is nothing but equihbrium 
condition of hydrostatic pressure: 

V(nor(x)) = (27) 

where no(r) — J dwf is the density determined from the solution. 

When we keep the constraint of the static solubihty conditions (16) and (18), there is no 

room to appear q*^^) in the hydrodynamic equations. Even when the perturbation violates 
the static solubility condition, the possible second order heat current is only 

in the perturbation, where o; = 45/8, and is the shear viscosity. It should be noted that the 
above expression is from eq.(15.3.6) in Chapman- Cowling [1] with neglecting the nonlinear 
terms of the velocity field. We also note that the perturbed pressure tensor has the form 
of Navier-Stokes order, because the velocity field itself is a perturbed quantity. Thus, the 
linearized hydrodynamics has a similar form as that of Navier-Stokes equation except for a 
term of the heat current in eq.(28). 

Nevertheless, it is still not easy to discuss the complete stability of hydrodynamic modes 
around the steady Burnett solution. Let us discuss the perturbation of density and temper- 
ature 

n = no(l + p(r, t)), O = T{x){l + ^(r, t)), (29) 

where is the temperature after we include the perturbation. Here we keep that the 
macroscopic velocity does not exist for the sake of simplificity. In this case it is easy to 
confirm that the linearized hydrodynamic equations are reduced to dtp — and 
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where the heat conductivity is assumed to be k = n* with K* = 756ii/(64(720F^). This 
equation apparently shows that the solution of 6 is relaxed to zero for all wavelength, as 
time goes on. Thus, the steady solution of Boltzmann equation at Burnett order is stable 
for the perturbation (29) without macroscopic flow. 

IV. MOLECULAR-DYNAMICS SIMULATION 

In this section, we show our preliminary result of two-dimensional molecular-dynamics 
(MD) simulation to examine the Knudsen effect under the heat conduction. Although we 
adopt two-dimensional systems, MD can check the validity of their implicit assumptions of 
the calculation by Kim and Hayakawa [10] where their calculation neglect the effects of the 
thermal wall at x = and use the distribution function at x = in eq. (23). 

Wc prepare two situations where Ti/Tq = 4 and Ti/Tq = 1/4. Number of total particles 
is 10,000 for both cases. Neglecting the non-Gaussian correction in velocity distribution 
function, the mean- free paths I — l/{2\/2na) for the center regions of four cells are respec- 
tively evaluated as Ineq = 11.24(/ and leg = 17.06d for the case Ti/Tq = 1/4 and Ineq = 21.40(i 
and /eg = 9.534(i for Ti/Tq = 4, where leq and Ineq represent the mean-free paths for an 
equilibrium cell and a nonequilbrium cell, respectively. Figure 2 shows the result of both 
situations. Lq ~ 1071d and the height H ~ 268. 6(^. The linear dimension along x-axis of the 
hole is 0.40845 d and the diameter of the hole is l.llSc?. This system has the area fraction 
= 0.02730. The vertical axis in Fig.2 plots the pressure divided by the density and tem- 
perature of equilibrium cell Tq. The horizontal axis shows the time of simulation measured 
by the number of collisions for each particle. We adopt the diffusive boundary condition for 
the walls at x — and x — ±Lo, and the periodic boundary condition at y = ±H/2. 
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FIG. 2. The time evolution of pressures in ]VID simulation where Pneqi 

and Peql are the result 
of the simulation for Ti/Tq = 4, and Pneq2 and Peq2 are the result for Ti/Tq = 1/4. 

In Fig. 2, the suffices 1 and 2 in the figure correspond to the cases of Ti/Tq = 4 and 
Ti/Tq = 1/4. Our result is contrast to those of both SST [9] and the calculation by Kim 
and Hayakawa [10]. Here, the osmosis AP depends on the gradient of temperature in which 
AP becomes positive for Ti/Tq = 4 and AP becomes negative for Ti/Tq = 1/4. This 
simulation suggests that both SST and the calculation by Kim and Hayakawa cannot be 
used for Knudsen effect in physical situations. 

There are two possibilities to have linear dependence of the osmosis AP on the tem- 
perature gradient: (i) Knudsen's layer effect is important on the thermal wall, and (ii) the 
assumption of small Knudsen number is violated. As a result of finiteness of Knudsen num- 
ber which is around 1/50-1/100, the temperature at the place deviated from the hole by 
the mean-free path is significantly different from the temperature at the hole. We believe 
that the second possibility plays a dominant role, because the gap of the pressure is en- 
hanced when the nonequilibrium cell has higher temperature, i.e. larger mean-free path 
Ineq = 21.40(i than that for low temperature nonequilibrium cell l^q = 11.24(i. The details 
of MD simulation will be reported elsewhere. 
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V. CONCLUSION 



In conclusion, we have explained our recent finding of a new explicit solution of Boltz- 
mann equation under the steady heat conduction. We apply the solution to examine the 
validity of SST and information theory, and find that both of them are not appropriate 
as they are. We also discuss the stabihty of the steady solution we have obtained, and 
confirmed stabihty at least for the case of no systematic fiow. Through MD simulation of 
Knudsen effect under the heat conduction, we have obtained the result that AP depends on 
the direction of the temperature gradient. This indicates that our naive application of the 
steady solution of Boltzmann equation to Knudsen effect is not accurate. This also suggests 
that the basic assumption of SST cannot be valid if the influence of physical wall exists. 
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TABLES 

TABLE I. The number constants 6or- in eq.(14) 



r 


bor 


blr 


&2r 





- 


- 


-3.320 X 10-2 


1 


- 


1.025 


-1.276 X 10-1 


2 


4.380 X 10"^ 


4.892 X lO"''^ 


6.414 X 10^2 


3 


-5.429 X 10-2 


3.715 X 10-3 


5.521 X 10-3 


4 


-4.098 X 10-^ 


2.922 X 10-^ 


4.214 X 10-3 


5 


-3.184 X 10""^ 


2.187 X 10-5 


3.106x10-5 


6 


-2.087 X 10-5 


1.492 X 10-6 


1.861 X 10-6 


7 




8.322 X 10-^ 
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